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A dynamical-system picture of a simple 
branching-process phase transition 

David Williams 



Abstract 

This paper proves certain results from the 'appetizer for non-linear 
Wiener-Hopf theory', 0). Like that paper, it considers only the simplest 
possible case in which the underlying Markov process is a two-state 
Markov chain. Key generating functions provide solutions of a simple 
two-dimensional dynamical system, and the main interest is in the way 
in which Probability Theory and ODE theory complement each other. 
No knowledge of either ODE theory or Wiener-Hopf theory is assumed. 
Theorem 1.1 describes one aspect of a phase transition which is more 
strikingly conveyed by Figures 14.11 and 14.21 

AMS subject classification (MSC2010) 60J80, 34A34 



1 Introduction 

This paper is a development of something I mentioned briefly in talks 
I gave at Bristol, when John Kingman was in the audience, and at the 
Waves conference in honour of John Toland at Bath. I thanked both 
John K and John T for splendid mathematics and for their wisdom and 
kindness. 

The main point of the paper is to prove Theorem 11.11 and related 
results in a way which emphasizes connections with a simple dynamical 
system. The phase transition between Figures 14.11 and 14.21 looks more 
dramatic than the famous 1-dimcnsional result we teach to all students. 
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The model studied here is a special case of the model introduced in 
. I called that paper, which contained no proofs, an 'appetizer'; but 
before writing a fuller version, I became caught up in Jonathan Warren's 
enthusiasm for the relevance of complex dynamical systems (in C 2 ). See , 
Warren and Williams This present paper, completely independent of 
the earlier appetizer and of my paper with Warren, can, I hope, provide 
a more tempting appetizer for what I called 'non-linear Wiener-Hopf 
theory'. No knowledge of any kind of Wiener-Hopf theory is assumed 
here. 

I hope that Simon Harris and I can throw further light on the models 
considered here, on the other models in Q, and on still other, quite 
different, models. 

Our model. A particle moving on the real line can either be of type + in 
which case it moves right at speed 1 or of type — in which case it moves 
left at speed 1. 

Let q- and q+ be fixed numbers with q_ > q + > 0, and let j3 be a 
positive parameter. We write K± = q± + f3. So, to display things, we 
have 

0_ > q+ > 0, > 0, K+ = q++ /?, K_ =q-+P. (1.1) 
We define 

/3c:=|(^-V9T) 2 . 

A particle of type ± can flip to the 'opposite' type at rate q± and can, 
at rate /3, die and at its death give birth to two daughter particles (of 
the same type and position as their 'parent'). This is why /? is a 'birth 
rate'. The usual independence conditions hold. 

Theorem 1.1 Suppose that our process starts at time with just 1 
particle of type + at position 0. 

(a) Suppose that /3 > f3 c . Then, with probability 1, each of infinitely 
many particles will spend time to the left o/O. 

(b) Suppose instead that /3 < j3 c . Then, with probability not less than 
1 — \/q+/Q—i there will never be any particles to the left o/O. 

Large-deviation theory (of which the only bit we need is proved here) 
allows one to prove easily that if (3 < f3 c , then, almost surely, only a 
finite number of particles are ever to the left of 0. 

The interplay between the Probability and the ODE theory is what is 
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most interesting. We shall see that j3 c plays the role of a critical para- 
meter in several ways, some probabilistic, some geometric. The 'balance' 
which occurs when (3 = f3 c is rather remarkable. 

The paper poses a tantalizing problem which I cannot yet solve. 

2 Wiener— Hopferization 

2.1 The processes {^(tp) : (p > 0} 

For any particle i alive at time t, we define to be its position on the 
real line at time t, and we extend the definition of <£>i by saying that at 
any time s before that particle's birth, 3>j(s) is the position of its unique 
ancestor alive at time s. 

So far, so sane! But we are now going to Wiener-Hopfcrize everything 
with a rather clumsy definition which defines for each if > two subsets, 
S + ((p,f3) and S~(<p,/3), of particles. 

We put particle i in set S + (ip,/3) if there is some t in \B(i),D(i)) 
where B(i) and D(i) are, respectively, the times of birth and death of 
particle i, such that 

• = <P, 

• ®i(t) > max{$j(s) : s < t}, and 

• i>i grows to the right of t in that, for e > 0, there exists a 6 with 
< 6 < e such that > ip throughout (t, t + 5). 

At the risk of labouring things, let me describe S~ (ip, (i) for <p > 0. 
We put particle i in set S~(ip,(3) if there is some t in [B(i), D(i)) such 
that 3>i(t) = —<p, &i(t) < min{$.j(s) : s < t}, and $ 2 ; decreases to the 
right of t in that, for e > 0, there exists a 6 with < 5 < e such that 
$,(•) < —ip throughout (t,t + 5). 

Of course, there may be particles not in \J V>0 {S + (ip, (3) U S~(ip 7 /?)}. 

We define N + (tp,f3) [respectively, N~(ip,j3)] to be the number of 
particles in 5 l+ (i y 9,/3) [resp., S~(ip,/3)]. 

We let Pg [respectively, P^] be the probability law of our model when 
it starts with 1 particle of type + [resp., — ] at position at time 0; and 
we let El" [resp., Er] be the associated expectation. 

We often suppress the '/3' in the notation for , E^ 1 , , N^ 1 . 

Then, under P + , N + = {N + ((p) : ip > 0} is a standard branching pro- 
cess, in which a particle dies at rate K + and is replaced at the '$-time' 
of its death by a random non-negative number, possibly 1 and possibly 



4 David Williams 

oo, of children, the numbers of children being independent, identically 
distributed random variables. I take this as intuitively obvious, and I am 
not going to ruin the paper by spelling out a proof. 

Note that in the P~ branching process N~ = {N~((p) : tp > 0}, a 
particle may die without giving birth. 

For < 6 < 1, dchnc 

g ++ (ip, 9) := E+9 N+ ^\ h~ + {y, 6) := E~9 N+ ^\ 
Clearly, for < 6 < 1, 

h- + (ip,e) = E-E-[^ + ^ | N+(0)] = E-g+ + (tp,e) N+ W 
= H-+(g++(<p t 0)), 

where 
where 

h~ + :=P"[A^ + (0) = n}. 
It may well be that := P"[iV+(0) = oo] > 0. Note that 
H- + (l-) := \\mH-+{9) = P~[iV + (0) < oo]. 



2.2 The dynamical system 

We now take 9 in (0, 1) and derive the backward differential equations 
for 

x^):=g++&,9), y(if):=h- + (<p,e), 

in the good old way in which we teach Applied Probability, and then 
study the equations. 0] looks a bit more 'rigorous' here. 

Consideration of what happens between times and dt tells us that 

x((p + dfp) = {1 - K+d<p}x(<p) +{q + dip}y(ip) + {P d<p}x(ip) 2 +o(d<p). 

The point here is of course that if we started with 2 particles in the + 
state, then E9 N ^ = x(ip) 2 . We see that, with x' meaning x'(ip), 

x = q + (y - x) + (i{x 2 - x). (2.1a) 

Similarly, remembering that $ starts to run backwards when the particle 
starts in state — , we find that 

y(ip - dip) = {1 - K- dtp}y(<p) + {q- d(p}x(ip) + {(3 d(p}y(<p) 2 + o(dip), 
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whence 

-y' = q-(x - y) + P{y 2 - y). (2.1b) 

Of course, y = H '"(x) must represent the track of an integral curve 
of the dynamical system (|2.1[) . and since y' = H~+ (x) x 1 , we have an 
autonomous equation for H h which we shall utilize below. 

Note that the symmetry of the situation shows that x = (y) must 
also represent the track of an integral curve of our dynamical system, 
though one traversed in the l ^p-TeveTsed , direction. 

Probability Theory guarantees the existence of the 'probabilistic solu- 
tions' of the dynamical system tracking curves y = H '"(a:) and x = 
H + -(y). 

Lemma 2.1 There can be no equilibria of our dynamical system in the 
interior of the unit square. 

Proof For if (x, y) is in the interior and 

q+{y -x) + P{x 2 -x) = 0, q-(x - y) + (3(y 2 - y) = 0, 

then y > x from the first equation and x > y from the second. Hence 
x = y and x 2 — x = y 2 — y = 0. □ 

2.3 Change of 

We need to think about how a change of would affect things. Suppose 
that a = E+0 Ar+ W') where < a < 1. Then 

E+E+[e N+ ^+^ I N+(<pj\ = E+a N+ ^ = g ++ (v,a), 

where a = g ++ {ij), 9). So, we have the probabilistic-flow relation 

g++(ip + iP,6) = g++(^g++^ 7 0)). (2.2) 

Likewise, h ^((p + tp,6) = h ^ (ip, g ++ (ijj, 6)). Thus, changing from 9 
to a = E + 9 N W just changes the starting-point of the motion along 
the probabilistic curve from (9,H ^(9)) to the point (a, H h (a)) still 
on the probabilistic curve. This is why we may sometimes seem not to 
care about 9, and why it is not in our notation for x(ip), y(<p). But we 
shall discuss 9 when necessary, and the extreme values and 1 of 9 in 
Subsection 14.61 

If for any starting point Vo = (xo,yo) within the unit square, we 
write V(< / 9, vq) for the value of (x(ip), y(<p)), then, for values of <p and ip 
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in which we are interested, we have (granted existence and uniqueness 
theorems) the ODE-flow relation 



which generalizes (|2.2|) . (The possibility of explosions need not concern 
us: we are interested only in what happens within the unit square.) For 
background on ODE flows, see [f|. 

3 How does ODE theory see the phase transition? 

3.1 An existence theorem 

Even if you skip the (actually quite interesting!) proof of the following 
theorem, do not skip the discussion of the result which makes up the 
next subsection. 

Theorem 3.1 There exist constants {a n : n > 0} with ao = 0, all 
other a n strictly positive, and 



and a solution (x(tp),y((p)) of the 'tp-reversed' dynamical system 



-x' = q+(y - x) + P{x 2 - x), y' = q_(x - y) + (3(y 2 - y), 

such that x(ip) = A(y(tp)) 7 where we now write A(y) = a n y n . 

Proof Assume that constants described exist. Since x'((f) = 

A '{y{f))y'(f), w e have 



- {q+y + (3A(y) 2 - K+A(y)} = A'{y) {q-A(y) + py 2 - K.y) . 



Comparing coefficients of y , 

/3oq - K + a = -aiq-aQ, 

and we are guaranteeing this by taking ao = 0. Comparing coefficients 
of y , we obtain 



V(yj + ip, v ) 



V(^,V(^v )) 




(3-1) 



q-a\ - {K- + K+)ai +q+ = 0. 



We take 



K- + K+ - y/(fr_ + K+Y - Iq~q~; 

2q- 
2q+ 



K- + K+ + sj (K- +K + ) 2 -4o_q + 
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from which it is obvious that < a\ < 1. 

On comparing coefficients of y n we find that, for n > 2, 

{K + + nK- — (n + l)<j_ai}a„ 

n— 1 n—2 

= j3 OAOn-ft + 9- + l)a fc+ ia n _ fc + /3(n - l)a„_i. 

fe=i fe=i 

We now consider the a„ as being defined by these recurrence relations 
(and the values of ao and ai). It is clear that the a„ are all positive. 
Temporarily fix N > 2, and define 

iV 

^iv(y) := 5Z a " y ™' 

n=0 
N 

L(y) := ^£„y™ := -g+j/ - /L4 N (y) 2 + 

:= 5> n y n := A' N (y){q_A N (y) + [3y 2 - K_y}. 

For n < N we have £„ = r n by the recurrence relations. It is clear that 
for n > N we have £ n < and r„ > 0. Hence, for all y in (0, 1), 

-q + y-/3A N (y) 2 +K + A N (y) < A' N (y){q_A N (y)+f3y 2 -K_ y }. (3.2) 

Suppose for the purpose of contradiction that there exists yo in (0,1) 
with A N (y ) = y . Then 

-q+ - fSy + K+ < A' N (y a ){q. + f3y a - K_}. 

However, the left-hand side is positive while the right-hand side is neg- 
ative. 

Because An(0) — and A' N (0) = a% < 1, the contradiction establishes 
that Ajy(y) < y for y G (0, 1), so that An(1) < 1. Since this is true for 
every N, and each a n (n > 1) is strictly positive, we have Ajv(l) < 1 for 
every N. 

By inequality (|3.2p . we have 

D N q_(A N - I) + q+ + 0A% - K+A N > 0, 

where An := Ajy(l) < 1 and Djv := A' N {1). Because each a n (n > 0) is 
positive it is clear that An < Dm- Wc therefore have 

q+ + f3A 2 N - K+A N > D N q_{l - A N ) > q-A N (l - A N ), 

which simplifies to 

{1 - A N ){q+ - {f3 + q_)A N } >0. 
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Since (1 — An) > 0, we have (ft + q^)A^ < q+, and result (|3.1[) follows. 
It is clear that we now need to consider the autonomous equation for 

y = y{<p)- 

y' = q-[A{y)-y]+f3[y 2 - y], y(0) = 6. 

But we can describe y(<p) as E6 Z ^ where {Z(tp) : tp > 0} is a classical 
branching process in which (with the usual independence properties) a 
particle dies at rate K_ and at the moment of its death gives birth to 
C children where 



P(C = n) = 



q-a n /K- if 1 < n < oo and n^2, 

(ft + q-a 2 )/K_ ifn = 2. 



Of course, aoo = 1 — a„. 

Then (x(<p),y(tp)) = (A(y(9?)), y((p)) describes the desired solution 
starting from (A(6),9). □ 



3.2 Important discussion 

Of course, ODE theory cannot see what we shall see later: namely that 
A(-) = #+-(•) when p > p c but A(-) ^ H+-(-) when /3 < /3 C . When 
/3 > /3 C , the curve a; = i? H (y) is the steep bold curve x = A{y) at 
the left-hand side of the picture as in Figure ETT1 But when ft < /3 C , the 
curve x = i?" 1 (y) is the steep bold curve at the right-hand side of the 
picture as in Figure H31 Ignore the shaded triangle for now. 

What ODE theory must see is that whereas there is only one integral 
curve linking the top and bottom of the unit square when ft > ft c , there 
are infinitely many such curves when ft < ft c of which two, the curves 
x = A(y) and x = (y), derive from probability generating functions 
(pgfs). 

It does not seem at all easy to prove by Analysis that, when ft < ft c , 
there is an integral curve linking the bottom of the unit square to the 
point (1, 1), of the form x = F(y) where F is the pgf of a random variable 
which can perhaps take the value oo. Methods such as that used to prove 
Theorem 13 . 1 1 will not work. 

Moreover, it is not easy to compute (0) when ft is equal to, or 
close to, ft c . If for example, q+ = 1, q_ = 4 and ft = 0.4, then one can 
be certain that (0) = 0.6182 to 4 places, and indeed one can easily 
calculate it to arbitrary accuracy. But the critical nature of ft c shows 
itself in unstable behaviour of some naive computer programs when ft is 
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equal to, or close to, j3 c . I believe that in the critical case when q + = 1, 
q- = 4 and /3 = 0.5, H + -(0) is just above 0.6290. 

Mathematica is understandably extremely cautious in regard to the 
non- linear dynamical system (|2.1[) . and drives one crazy with warnings. 
If forced to produce pictures, it can produce some rather crazy ones, 
though usually, but not absolutely always, under protest. Its pictures 
can be coaxed to agree with those in the earlier appetizer which were 
produced from my own 'C Runge-Kutta program which yielded Post- 
script output. Sadly, that program and lots of others were lost in a 
computer burn-out before I backed them up. 

4 Proof of Theorem 11.11 and more 

4.1 When (3 > (3 C 
Lemma 4.1 When /? > /3 C , 

(a) H-+{0) = p-[JV+(0) = 0] = 0, 

(b) H-+(l-) = F-[N+(0) < oo] < 1, 

(c) H+-(l-) = P+[7V-(0) < oo] < 1, 

(d) H+-(0) = P+[JV-(0) = 0] = 0. 

It is clearly enough to prove the lemma under the assumption 

\{^-s/q+) 2 <P< i(v^T+\/9+) 2 , 

and this is made throughout the proof. 

Proof Result (a) is obvious. 

The point (1,1) is an equilibrium point of our dynamical system, and 
we consider the linearization of the system near this equilibrium. We 
put x = 1 + £, y = 1 + ?/ and linearize by ignoring terms in £ 2 and rf: 

-q+ + P q+ \ ft\ 
-q- q- - PJ \vj ' 

the matrix being the linearization matrix of our system at (1,1). The 
characteristic equation for the eigenvalues of this matrix is 

A 2 + (q+ - q_)\ + (q_ + q + )(3 - /3 2 = 0. 

The discriminant l B 2 — \AC is 

{2[3-{q + +q-)} 2 -Aq_q + . 
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This expression is zero if /3 = \{^Jq~^ ± \/9+) 2 ' m our case i the 
eigenvalues A have non-zero imaginary parts. Any solution of our system 
converging to (1, 1) as ip ±oo must spiral, and cannot remain inside 
the unit square. Hence results (b) and (c) hold. 

It is now topologically obvious (since there are no equilibria inside the 
unit square) that we must have P + [N~ (0) = 0] = 0; otherwise how could 
the curve x = (y) link the top and bottom edges of the square? Thus 
result (d) holds. □ 

Of course, we can now deduce from result (|3.ip that (when /3 > /3 C ) 

H+-{1-) = P+[JV-(0) < oo] < q+/(q- + 13). 

Figure 14.11 which required 'cooking' beyond choosing different ip- 
rangcs for different curves, represents the case when q + = 1, g_ = 4 



Dynamical- system picture of a phase transition 



11 



and j3 = 4. The lower bold curve represents y = H h (cc) and the upper 
x = (y). As mentioned previously, (y) = A(y) where A is the 
function of Theorem 13.11 

The motion along the lower probabilistic curve y = H '"(x) will start 
at (8,H and move towards (0,0) converging to (0,0) as ip —> oo 

since N ((f) — > oo. If wc fix ip and let 9 1, we move along the curve 
towards the point (1,H(1— )). Of course, we could alternatively leave 
fixed and run tp backwards. It is clear because of the spiralling around 
(1,1) that the power series H ^(x) must have a singularity at some 
point x not far to the right of 1. 

Motion of the dynamical system along the steep probabilistic curve 
x = (y) on the left of the picture will be upwards because it is the 
(^-reversal of the natural probabilistic motion. Now you understand the 
sweep of the curves in the top-right of the picture. 



4.2 A simple large-deviation result 

Let {X(t) : t > 0} be a Markov chain on {+, — } with Q-matrix 

-q+ q+ 
q- -q- 

Let V be the function on { + , — } with V(+) = 1 and V(— ) = —1 and 
define &x(t) = I V(X s )ds. Almost surely, <&x(t) -> oo. We stay in 
'dynamical-system mode' to obtain the appropriate Feynman-Kac for- 
mula. 

Let fi > 0, and define (with the obvious meanings of K ± ) 

u(t) = E+ cxp{-/i$ x (t)}, v{t) = E~ cxp{-^x(t)}- 

Then 

u(t + dt) = {l-g + dt}e _Mdt u(t) + g + dtu(t)+o(dt), 
with a similar equation for v. We find that 
u \ f^q+^t 1 q+ \( u \ ( u {t)\ r irt \r\~\ 



v' J V q_ — g_ + /jy \v J ' \v(t) 



= exp{t(Q- M y)} I 



where V also denotes the operator ( Q _ 1 ) of multiplication by the func- 
tion V. 

Lemma 4.2 If fi < j3 c :— \ (jqZ— \fq+) 2 , then there exist positive 
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constants e, k, A such that 

e /3 *P ± [$x(i) < et] < Ae~ Kt . 
Proof We have just shown that 

E ±[e-M**« /(Xt) ] = exp{t(Q - fiV)}f. 
Now Q — \iV has larger eigenvalue 



7 = -W- + q+) + W(l- + 1+) 2 - 4 (<?- - 9+) A» + V- 

We fix /x at g(<7- — <?+) to obtain the minimum value | (^/gT — v / 9+) 2 
of 7. Hence, for e > and some constant ^4 £ , 

P ± [$x(f)<rf] = P ± [n(et- $x(t)) > 0] < E ± exp{/iei- 

< ^ e exp|i£(g_ - - \ (y/ql- V«+) 2< } • 

The lemma follows. □ 

For a fine paper proving very precise large-deviation results for Markov 
chains via explicit calculation, see Brydges, van der Hofstad, and Konig 

!• 

4.3 When p < p c 
Lemma 4.3 When ft < f3 c , 

(a) H-+(0) = P-[2V+(0) = 0] = 0, 

(b) H-+(l-) = p-[JV+(0) < 00] = 1, 

(c) H+-(0) = P+[JV-(0) = 0] > 0, 

(d) H+-(l-) = P+[AT-(0) < 00] = 1. 

Note Though Figure B~2l relates to the case when (3 = f3 c , pictures for 
the subcritical case when j3 < j3 c look very much the same. 

Proof Result (a) remains trivial. 

By Lemma [4.21 there exist e > 0, n > and A > such that, for a 
single particle moving according to Q-matrix Q, we have 

e pt ¥>[$x(t) < et] < Ac~ Kt . 

For the branching process, the expression on the left-hand side is the 
expected number of particles with $-value less than or equal to et at 
real time t. So the probability that some particle has $-value less than 
or equal to t is at most Ae~ Kt . 
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Figure 4.2 Mathematica picture of the critical case when q+ — 1, 
q- = 4, 13 = 0.5 

By the Borel-Cantelli Lemma, there will almost surely be a random 
positive integer no such that for all n > no, every particle alive at real 
time n will have $- value greater than en. Since $ can only move left at 
speed 1, there must almost surely come a time after which no particle 
has a positive $- value. Hence P~[.ZV + (0) = oo] = 0, and result (b) is 
proved. 

Now suppose for the purpose of contradiction that P + [iV~ (0) > 0] = 1. 
Since a particle started at state + can remain there for an arbitrary long 
time without giving birth, it follows that any particle in the + state and 
with any positive $- value will have a descendant for which $ will become 
negative. This contradicts what we proved in the previous paragraph, so 
result (c) is established. 

Since the y = H '"(x) curve connects (1,1) to (0,0) and the other 
probabilistic curve x = (y) starts at (7J" 1 (0),0) where (0) > 
0, and since these curves cannot cross at an interior point of the unit 
square, it must be the case that H A (1 — ) = 1, so that property (d) 
holds. □ 
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In the analogue of Figure 14.21 for a sub-critical case (which, as I have 
said, looks very much like Figure FT~2"j) . motion along the higher prob- 
abilistic curve y = H h (x) will again start at (6, H ^(0)) and move 
towards (0,0), this because N + (ip) — > oo. Since N~(tp) — > 0, the natural 
probabilistic motion of the lower curve x = (y) will converge to 
(1,1); but the tp- reversal means that the dynamical system will move 
downwards along this curve. 

Sketch of geometric proof that H h (l — ) = 1 if (3 < (3 C It is enough to 
prove the result when (3 = j3 c . Let m = \/q-/q+, the slope of the unique 
eigenvector of the linearity matrix at (1, 1). Draw the line of slope m 
from (1, 1) down to the y-axis, the sloping side of the shaded triangle 
in the picture. Now it is particularly easy to check that at any point of 
the sloping side the (dy/da;)-slope of an integral curve is greater than 
m. If the convex curve y = H ^(x) hit the vertical side of the triangle 
at any point lower than 1, we would have 'contradiction of slopes' where 
it crossed that sloping side. □ 

4.4 Nested models and continuity at phase transition 

Take /?o > f3 c and let Mp be our model with initial law Ft (in the 
obvious sense). Label birth-times T\, T2, T3, . . . in the order in which 
they occur, and for each n call one of the two children born at T n 'first', 
the other 'second'. Let U\, U2, U$, ... be independent random variables 
each with the uniform distribution on [0, 1). We construct a nested family 
of models {Mp : (3 < (3q} as follows. 

Fix (3 < /3q for the moment. If U„ > erase the whole family 

tree in Mp B descended from the second child born at time T n . Of course, 
this family tree may already have been erased earlier. In this way, we 
have a model Mp with desired law P^~. The set S~{(p 1 (3) will now denote 
the S~(<p) set for the 'nested' model Mp, and N~(ip,(3) will denote its 
cardinality. 

A particle i contributing to S~(Q,/3) determines a path in {4-,—} x 
[0,oo): 

{ (Ancestor, (t), $ 4 (t)) : t < p t } 

where pi is the first time after which $^ becomes negative. Along that 
M^-path, there will be finitely many births. Now, for fixed (3 it is almost 
surely true that U n 7^ j3 for all n. It is therefore clear that, almost surely, 
for f}' < (3 and f3' sufficiently close to /3, the M^-path will also be a path 
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of Mpi. In other words, we have the left-continuity property 
S ,_ (0,/3)= (J S-(O,0'), almost surely. 

i8'<0 

It therefore follows from the Monotone-Convergence Theorem that 

E + N~(0,/3 C ) =f UmE + JV _ (0,/3). (4.1) 

Clearly, something goes seriously wrong in regard to right-continuity 
at /3 C . Suppose we have a path which contributes to S~(0,/3) for all 
f3 > f3 c . Then, for all birth-times T n along that path we have U n < 
for all /J > /3 C and hence U n < /3 c //?o- Hence 

S-(p,p c )= f) S"(0,/3). 

/3>/3 c 

But it is clearly possible to have a decreasing sequence of infinite sets 
with finite intersection. And recall that (more generally) the Monotone- 
Convergence Theorem is guaranteed to work 'downwards' (via the Dom- 
inated- Convergence Theorem) only when one of the random variables 
has finite expectation. 

4.5 Expectations and an embedded discrete-parameter 

branching process 

If either of the curves y = H h (x) or x = (y) approaches (1,1), 
it must do so in a definite direction and it is well known (and an im- 
mediate consequence of l'Hopital's rule) that that direction must be an 
eigenvector of the linearity matrix at (1,1). When (3 = (3 C , there is (as 
we have seen before) only one eigenvector ( ™ ) with m = \/q-/q+- Thus 

E-N+(0, &) = (H-+)'(1,/3 C ) =m= y/q-/q+. (4.2) 

We know that if f3 < /? c , then (1— ) = 1 and we can easily check 
that (as geometry would lead us to guess) 

E+N-(0,/3) = (H + -)'(1,P) < l/m = y/q+/q-, 

and now, by equation (|4.ip we see that 

E+N- (0, (3 C ) < l/m = y/q+/q-. (4.3) 

In particular, P + [iV - (0, (3 C ) = oo] = 0, and so, in fact, (1) = 1 and 
we have equality at (|4.3p . whence 

P+[AT-(0,/3 c ) > 1] < E+AT-(0,/3 c ) = y/q-/q+, (4.4) 
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part of Theorem 11.11 
Now, for if > 0, let 



K<p) 



E + N + {lp) 7 c(ip) = E-N + {ip). 



Then 



q+Mp) - Hv)} + /?%), etc., 



so that the linearization matrix at (1, 1) controls expectations. We easily 
deduce the following theorem. 

Theorem 4.4 When = fi c , 



E+N-(ip) = fs± e -^ q -- g +^, E-N-(<p) = e~^ q -- q +^. 



For any /3, we can define the discrete-parameter branching-processes 
{W /± (n) : n > 0} as follows. Let Bi be the birth time and Di the death 
time of particle i. Recall that $j is defined on [0, Di). Let (T,(0) = and, 
for n > 1, define 

(7,(71) := inf{f :Bi<t<Di:t> a^n - 1); (-l)"$»(i) > 0}, 

with the usual convention that the infimum of the empty set is oo. Let 

W + {n) := : o;(2n) < oo}, W~(n) := : crj(2n + 1) < oo}. 

The 'W notation is suggested by 'winding operators' in linear Wiener- 
Hopf theory. 

Theorem 4.5 W is a classical branching process under Wt, and is 
critical when j3 = j3 c . 

The proof (left to the reader) obviously hinges on the case when tp = 
of Theorem 14.41 And do have a think about the consequences of the 
'balance' 





E-N+{ V7 f3 c )E + N-{v 7 (3 c ) = 1 



in that theorem. 



4.6 When 9 = or 1 



If we take 9 = and set 



K<p) ■= P + [N~ (tp) = 0] , c(<p) := P- [N- (tp) = 0] , 
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then {(b(tp) , c(ip)) : if > 0} is a solution of the ^-reversed dynamical 
system such that 

b(ip)=H+-(c(ip)). 

When f3 > f3 c , this solution stays at equilibrium point (0,0). When 
(3 < P c , {b(<p),c(<p)) moves (as < ip f oo) from (ff+-(0),0) to (1,1) 
tracing out the right-hand bold curve in the appropriate version of Figure 

m 

When 6=1, {(B(tp),C(tp)) : <p > 0}, where 

:= P+[iV-(^) < oo], C(<p) := ¥-[N~{ip) < oo], 

gives a solution of the (^-reversed dynamical system. When j3 < f3 c , this 
solution stays at the equilibrium point (1,1). When j3 > f3 c , (B(ip), C(ip)) 
moves (as < ip t oo) from (H^ (1— ), 1) to (0, 0) tracing out the bold 
upper curve in the appropriate version of Figure 14.11 Of course, there 
is an appropriate version ('with + and — interchanged') for the lower 
curve. 

4.7 A tantalizing question 

When (3 > f3 c , we have for the function A(-) = A(-,/3) of Theorem 13. II 

When P < f3 c , we have, for < 9 < 1, A(9,f3) = E+6 Y « for some 
random variable Yp. Can we find such a Yp which is naturally related to 
our model? In particular, can we do this when f3 = /3 C ? It would be very 
illuminating if we could. 

What is true for all (3 > is that if X and $x are as at the start of 
Subsection l4~2l and r^(0) := inf{t : $ x (t) < 0}, then 

E+ expj-^O)} = oi, 

with a± as in Theorem 13.11 and this tallies with a\ = P + [A r_ (0) = 1] 
when (3 > f3 c . Proof of the statements in this paragraph is left as an 
exercise. 
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